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. Abstract 

, The fully compressible semi-geostrophic system is widely used in the modelling of 

large-scale atmospheric flows. In this paper, we prove rigorously the existence of weak 
Lagrangian solutions of this system, formulated in the original physical coordinates. In 
addition, we provide an alternative proof of the earlier result on the existence of weak 
solutions of this system expressed in the so-called geostrophic, or dual, coordinates. 
The proofs are based on the optimal transport formulation of the problem and on 
• recent general results concerning transport problems posed in the Wasserstein space 

^ ' of probability measures. 

1 Introduction 
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^ ' The behaviour of the atmosphere is governed by the compressible Navier-Stokes equations, 

CO , together with the laws of thermodynamics, equations describing phase changes, source 

terms and boundary fluxes. These equations are too complex to be solved accurately in a 
large-scale atmospheric context and therefore reductions and approximations of the Navier- 
Stokes equations are often used to validate and understand the solutions that have been 
computed. 

One such approximation, valid on scales where the effects of rotation dominate the 
flow, is the system of semi-geostrophic equations, see equations p.ip - p.Sp below. First 
introduced by Eliassen [S] and then rediscovered by Hoskins [TT], the semi-geostrophic 
equations are an extremely useful model, particularly in describing the formation of fronts, 
' and are both studied theoretically and used widely for numerical modelling and simulations. 

^ , In this paper, we prove rigorously the existence of weak Lagrangian solutions of the 

fully compressible semi-geostrophic system. The proof combines ideas of several previous 
papers, starting with the pioneering work of Benamou and Brenier [3] on the formulation of 
the semi-geostrophic equations as an optimal transport problem, and modern methods in 
the analysis of Hamiltonian ODEs in probability spaces, in particular the work of Ambrosio 

For an accurate representation of the behaviour of large-scale atmospheric flow, one 
should consider the fully compressible semi-geostrophic equations with variable Coriolis 
parameter and a free surface condition. The complexity of this problem means that so far 
results are obtained after relaxing one or more of these conditions. 

In [4] , Benamou and Brenier assumed the fluid to be incompressible, the Coriolis param- 
eter / constant and the boundaries rigid. They then used a change of variables, introduced 
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by Hoskins in [TT], to derive the so-called dual formulation and to prove that stable weak 
solutions of the semi-geostrophic equations in these dual variables exist, by interpreting the 
equations as a coupled Monge- Ampere transport problem. 

In [7] , Cullen and Gangbo relaxed the assumption of rigid boundaries with a more phys- 
ically appropriate free boundary condition. However, they additionally assumed constant 
potential temperature to obtain the 2-D system known as the semi-geostrophic shallow wa- 
ter system. After passing to dual variables, they showed existence of stable weak solutions 
for this system of equations. 

In [8 , Cullen and Maroofi proved that stable weak solutions of the 3-D compressible 
semi-geostrophic system in its dual formulation exist, returning to the assumption of rigid 
boundaries. 

The main problem posed by the existence results in [?], [7] and [5] is that they are all 
proved in dual space. It is difficult to relate these results directly to the Navier-Stokes 
equations, or indeed other reductions of them. For this reason, and in order to give the 
dual space results physical meaning, Cullen and Feldman [6^ mapped the solutions back to 
the original, physical coordinates and extended the results of [4] and [7], proving existence 
of weak Lagrangian solutions in physical space, in the incompressile case. 

In this paper, we utilise recent results in the analysis of ODEs in spaces of measures, 
in particular those of [5], [3], in order to provide an alternative proof of the dual space 
result in [5]. We also extend considerably the results of [S] to prove the existence of weak 
Lagrangian solutions of the fully compressible semi-geostrophic system in the original phys- 
ical coordinates. This is based on the existence of an appropriate flow map with rather low 
regularity; however we also show that, if we could assume additional regularity, then the 
solutions derived would determine classical solutions. 

The paper is organised as follows: in Section [5] we introduce various definitions, and the 
notation to be used throughout; in Section [3] we summarise the existing results regarding 
existence of solutions in dual space that are necessary for the theory presented in the 
following sections; in Section 2] we use the results of to provide an alternative proof 
of the existence of solutions in dual space; in Section [5] we deflne the concept of a weak 
Lagrangian solution and formulate our main theorem; in Section [S] we prove the existence 
and important properties of a Lagrangian flow in dual space; finally, in Section [7] we map 
the dual space Lagrangian fiow to physical space and use it to prove our main result, namely 
the existence of a weak Lagrangian solution. 

2 Useful Conventions, Notation and Definitions 

Physical quantities and constants (21) 

(i) Q, denotes an open bounded convex set in R'^, representing the physical domain con- 
taining the fluid; r > is a flxed positive constant; all functions in physical coordinates 
are deflned for {t, x) G [0, r) x il; 

(ii) u{t,x.) represents the 3-D velocity of the fluid; 

(in) u^(t,x) — (uf (t, x), x), 0) represents the geostrophic velocity; 

(iv) p{t,x.) represents the pressure; 

(v) /9(i,x) represents the density; 
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(vi) 6{t, x) represents the potential temperature. Given its physical meaning, we assume 
0{t,x) to be strictly positive and bounded; 

{vii) (/)(x) is the given geopotential representing gravitation and centrifugal forces. We 
assume that cj) G C^(r2) and that g|^(/)(x) ^ for all x e f2; 

(via) /cor denotes the Coriolis parameter, which we assume to be constant; 

(ix) pref is the reference value of the pressure; 

(x) Cy is a constant representing the specific heat at a constant pressure; 
{xi) Cp is a constant representing the specific heat at a constant volume; 

{xii) R represents the gas constant and satisfies R = Cp — Cy-, 

{xiii) K = Cp/cy denotes the ratio of specific heats (this is approximately 1.4 for air). 

Notations and other conventions (2-2) 

• Throughout, we will only consider measures that are absolutely continuous with re- 
spect to Lebesgue measure. Given an open set ^ in M^, we will denote by 

Pac{A) - the set of probability measures in with supports contained in A; 

Xa - the characteristic function of A. 

• Unless otherwise specified, measurable means Lebesgue measurable and a.e. means 
Lebesgue-a.e. 

• Dt denotes the Lagrangian derivative, defined as Dt = dt + u ■ S/. 

• The unit vector (0,0, 1) in is denoted by es. 

• For convenience, we will sometimes use the notation -F(t)(') = F{t,-) to denote the 
map F evaluated at fixed time t. In fact, we will often leave out time t as an argument. 

Important definitions 

Definition 2.1. We define 

P^,{R') := {H e Pac{R') ■■ I |x|V(x)dx < +(^}, 

with tangent space 

T^PliM.') = {V^:veC^{R^)f^'''^'\ (2.3) 

Definition 2.2. Let j be a probability measure on the product space x M^. Then, Borel 
probability densities jjL\ and 1x2 in R? are called marginals of 7 if 

Ml [A] =j[AxM.% fi2 [B] = 7[R3 X B] 

for all measurable A,BcM.^. We define 

r(/ii,/^2) = {7:7 probability measure on x with marginals /ii and ^2]- 
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Definition 2.3. Given two Borel probability densities /ii(-) and /i2(') in , we define the 
Wasserstcin-2 distance, W2, between fj,i and fi2 o,s follows: 

W^2(Mi,M2):= inf / |x-y|2d7(x,y). (2.4) 
We denote by Tq{pi, (12) the set of minimis ers of ^2^. 

The Wasserstein distance indeed defines a distance in the space of probability measures on 
M.^ (or more generally, any complete separable metric space); it can be used as an optimal 
transport cost between these measures, see |14) . 

Definition 2.4. Let F :MP ^M!^ be a measurable map. Let /ii be a measure in MP . Then 
we say that the measure ^2 on M"^ is the push forward of /ii with respect to F, denoted 
= fi2, if 

A*2[A] =^i[F-i(A)] 

for all measurable A C M"^. If, in addition, ni G L"'^(R'^), /i2 G L"'(M'^) and F^jii = fi2, 
then (see 16, Corollary A. 3]) 

[ (/;(F(x))/ii(x)dx= [ ip{X)^l2iX)dX, \f ifi e L°^(R^). (2.5) 

Remark 2.1. Since we are only interested in measures that are absolutely continuous with 
respect to Lebesgue measure, we can guarantee existence of a unique minimiser 70 in (12.4^ : 
see, for example, [3l Section 6.2.3]. We also have that 70 = (id, R^j)#/^i for some map 
R^J : R'^ — ?> which coincides /ii — a.e. with the gradient of a convex function. Thus, the 
map R^J is the unique minimiser of 

R^ / |x- R(x)|^ ^i(x)dx (2.6) 

JR3 

over all R such that R#/ii = fJ,2- 

Definition 2.5. Given a set A CZ MP we say that a sequence (nn) C PadA) is narrowly 
convergent to fi d Pac{A) as n ^ 00 if 

lim / (p{x) ^n{x) dx = I (p{x)ij{x)dx 
J A J A 

for every function ip £ Cb{A), the space of continuous and bounded real functions defined 
on A. 

Definition 2.6. Given a set A C M.^ and a function f : A ^ M, we say that f is semi- 
concave if there exists a constant k > such that the map x — > f{x) — ^k'^\x\'^ is concave in 
A. We call k a constant of semi-concavity. This definition can be generalised by replacing 
— jxp with any concave function. 

Note that if g is semi-concave, then Vg has bounded variation. 
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3 The semi-geostrophic equations: formulation and ex- 
isting results 

The fully compressible semi-geostrophic equations posed in the domain [0, r) x D,, with rigid 
boundaries dil, are the following system of equations (see, for example, 0): 

+ /cores X u + V0 + = 0, (3.1) 
P 

Dt9 = 0, (3.2) 

A- = -V • u, (t x) e [0, T)xn (3.3) 
P P 

/cores xu9 + V(/.+ -Vp = 0, (3.4) 
P 

P = Rp0(—] ' ; (3.5) 

u-n = 0, (t,x) G [0,t) X arj. (3.6) 

The unknowns in the above equations are = (uf,uf,0), u = (ui, 1*2,^3), p, p, 9. 

Equation p.ip is the momentum equation; p.2p represents the adiabatic assumption; 
(|3.3p is the continuity equation and p.4p represents hydrostatic and geostrophic balance. 
The equation p.Sp is the equation of state which relates the thermodynamic quantities to 
each other, and (j3.6p is the rigid boundary condition, where n is the outward normal to 

dn. 

The semi-geostrophic equations are a valid approximation when j^^^ << 1, and are 

accurate when -j- < where U is a typical scale for horizontal speed; L is a typical 

horizontal scale; iJ is a typical vertical scale; N is the buoyancy frequency. 

The energy associated with the flow, known as the geostrophic energy, is defined as 



i|u3|'(t,x)-f-0(x)+c,0(t,x) 

2 V Prof 



pit,x)dx. (3.7) 



3.1 Dual formulation 

In [5], solutions were obtained using a transformation into the so-called geostrophic, or 
dual, coordinates y = (yi , y2 , 2/3 ) • The coordinate transformation is given by: 

T : 17 ^ A C M^ T(t, x) = {T,{t, x), T2{t, x), T^it, x)) = (yi, 2/2, ys), 

with 

2/1 = a;i + /coJui(t,x), y2 = 3^2 - /cot"i(*> x), j/3 = 6'(i,x). (3.8) 

Note that, by ^^(vi), A C x [S, i] for some < S < I. 

Using this transformation, as well as (j3.5l) . it was shown in that we can write the 
energy in (|3.7I) as 

Ejt) ^ [ ^^-^'-^ - ^^'^ + - + ajt, X) dx K, fiait,x)rdx 
Jn 2/3 Jn 

f yiAl^, ri(t,x)p + |X2 - r2(t,x)p} + 0(x) /■ 

= / rr. N f7(t,x)dx + ifl / (Cr(t,x)) dx 

A2 T3(t,x) 

(3.9) 
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(\ K — 1 
I is constant and 
Prof/ 

a{t,x) ■.^0{t,x)p{t,x). (3.10) 

Using the results of 12^, it was shown that stable solutions of (I3.1l) - (j3.6p correspond to 
solutions that, at each fixed time t, minimise the energy E given by p.9p . We refer to this 
requirement as Cullen's stability condition. 

This stability condition can be formulated in terms of optimal transport concepts. In- 
deed, define the potential density v := T#(T as the push forward of the measure cr under the 
map T (see Definition l2.4l) . Then Cullen's stability condition is equivalent to the condition 
that, at each fixed time, the pair {a, T) minimises the energy p. 91) amongst all pairs (/i, T) 
with /i e Pac{^) and Tt^/i = h'. It was shown that, for a stable solution, the energy in (|3.9p 
can be written, dropping the dependence on time, in the following form 

E^{<7) ^ £{a, v)+Ki [ {ait, x))" dx, (3.11) 



where 

£{a,iy)^_M / c(x,T(t,x))cr(t,x)dx (3.12) 
with c a cost function, given by 



^fL{\^i-yi\' + \x2^y2?} + cl^{^) 



c(x,y)^ ^ — ' ■ - ^ -V (3^^3^ 

ys 

The energy minimisation required by Cullen's stability condition has now been refor- 
mulated as an optimal transport problem. Cullen's stability condition amounts to the 
requirement that the change of variables T{t, ■) from physical coordinates to geostrophic 
coordinates given by p.8p is the optimal transport map between cr ans i'. 

In 'W , it is shown that, given a potential density v, there exists a unique minimiscr a of 
p. lip . Given this pair ly, and cr, there always exists a unique minimiser T in p.l2p . given 
by the optimal map in the transport of a to v. the map T admits a unique inverse T~^, 
which is the optimal map in the transport of to cr with cost c(y, x) = c(x, y). 

Using the stability condition and the transformation given by p.Sp . the compressible 
semi-geostrophic equations p.ip - p.4l) can be written in dual variables as follows [S]: 

9tJ^(i,y) + V-(i^(i,y)w(t,y)) =0 in [0, r) x A, (3.14) 

w(t,y) =u^'(t,S(<,y)) = /cores x [y-S(i,y)], (3.15) 

S{t,-)=T-\t,-), (3.16) 

T(t, •) is the unique optimal transport map in p.l2p . (3-17) 

cr{t, •) minimises Ei,(^t,-)i') over Pad^), (3.18) 

jy(0, ■) = t/o(-) e ^'■(Ao), r e (1, oo), Ao C compact. (3.19) 



Equation p.l4p is the continuity equation satisfied by the potential density; p.isp 
defines the geostrophic velocity in dual variables; p.l6p . (I3.17P and p.lSp are required for 
a stable solution; and p.l9p is the prescribed initial condition. 



6 



The main result of [5] is the existence of a stable weak solution of the semi-geostrophic 
system in dual variables p.l4p - p.l9p . proved by an approximation procedure. In the same 
paper, several important results concerning optimal transport are proved. In Section [321 
we summarise these results. 

Remark 3.1. Since vq has compact support, for all t € (0,t) the solution J^(t,y) of the 
evolution p. 141) starting at time t = from i/q has compact support in R'^. By l\2.1^ (vi). 
supp{v) is contained in a bounded open set A, which is dependent on r, such that A C 
X [5, for some 5 with < (5 < 1. This follows from a standard fixed-point argument; 
see, for example, [SJ Section 3.5.2]. 

Remark 3.2. In [TU], Faria et al. extend the results of 6 to the case when the initial 
potential density vq is in . The question of whether this can be done for the fully 
compressible semi-geostrophic equations remains open. 

3.2 Existence and uniqueness for the energy minimisation 

Assume that 

• 51 and A are bounded open domains in R'^, 

• is convex ; A C x 

where the assumption on the vertical coordinate of A is justified by \2.1\i (vi). CuUen and 
Maroofi proved in [SJ Theorem 4.1] that, given v e Pac(A), there exists a unique a e Pac{^) 
that minimises the energy E^{-): 

Theorem 3.3. Assume that £7 and A satisfy I13.20\) . Let v e Pac(A) and assume that cj) 
satisfies V2. l]) (\yi). Then there exists a unique minimiser a of E^{-), over Pac{^)- This 
minimiser is given by iS.ll]) . 

By fixing this minimising cr, we can complete the energy minimisation by considering 
p.l2p as an optimal transport problem. Indeed, given two probability densities a G Pac{^) 
and v G Pac(A), consider the following optimal transport problem: 

_inf /^[T], with /<,[T] := / c(x,T(x))cr(x)dx, (3.21) 
T#T=i' Jn 



1 (3.20) 

S, - for some 6 with < S < 1, 

' 



where c(x, y), given by p.l3p . is the cost of transporting one unit of mass from x to y. 
The Kantorovich relaxation of p. 211) amounts to finding 7 G r(cr, v) that minimises 

/[7] := / c(x,y)d7(x,y), (3.22) 

with r(cr, v) given in Definition l2.2l Note that every transport map T generates a transport 
plan 77p G r((T, v) defined by 

7^:= (id,T)#a, 

such that /<,[T] =I[7^]. 

Now consider the Kantorovich dual problem corresponding to p.2ip (see, for example, 
[H Chapter 5]): 

sup J(a.^)if,g), with J(^_i,)(/,5) := / /(x)cr(x)dx+ / 5(y)i^(y)dy, (3.23) 
{f,g)eLip, Jn Ja 
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where 

Lip, : / £ W^^^{^1), g £ W^^^{K), /(x)+.g(y) < c(x,y) for all (x, y) G fix A}. 

The solution of the dual problem is a crucial ingredient of proving the main result of this 
section, Theorem 13.51 below . 

Let us discuss the dual formulation. We start by defining two operators that are standard 
in the theory of optimal transport. 

Definition 3.1. Let Q, and A be hounded open sets in R"^, and let c(x, y) he given hy \S.13\) . 
Assume that f is a function from Q into R and g is a function from A into R. Then the 
c-transform of f on is defined as 

ny) inl{c(x, y) - /(x)}, y e R-^ \ {y : yg = 0}, (3.24) 
xer2 

and the c-transform of g on A is defined as 

5e(x) inyc(x, y) - ^(y)}, x G R^. (3.25) 

yeA 

Cullen and Maroofi prove and utilise various useful properties of these c-transforms to 
show the existence of solutions of the Kantorovich dual problem. We summarise these 
results in the following (see [SI Lemma 3.1, Theorem 3.2] for proof): 

Theorem 3.4. Assume il and A satisfy i3.20]) . Let a E Pad^) and v G Pac(A). Assume 
that c(x, y) is given hy i3.13\) and that (f> satisfies i2.1\} fbi). Then the prohlem i3. 23\) of 
maximising J[cr,v) over Lip, has a solution {fo,go) such that {foY — go, {go)c = /o ^.i^d 

ll/olUi,=o(a) + ll5olUi,oo(A) ^C = Cin,A,c{; •)); 

(a) for V ~ a.e. y G A, the infimum in {3.24^ is attained at a unique point x = x(y) with 

V5o(y) = Vyc(x,y), (3.26) 

and for a — a.e. x G il, the infimum in LS. 25]) is attained at a unique point y — y(x) 
with 

V/o(x)- Vxc(x,y); (3.27) 

(Hi) go is semi-concave in A, /g is semi-concave in il and the constant of semi-concavity 
in hoth cases depends only on fl, A and c(-, •). 

Finally, we state the main existence result for the optimal transport problem, see [5J 
Theorem 3.3] for a proof. 

Theorem 3.5. Assume that ft and A satisfy \3.20\) . Let a G PaciS^ and v G Pac{A). 
Assume that 4> satisfies \2.1\) (i'yi). Then there exist maps T : fi A and S : A — 51, 
unique a — a.e. and v — a.e. respectively, such that 

(i) T is optimal in the transport of a to v with cost c(x, y), 

(ii) S is optimal in the transport of v to a with cost c(y,x) = c(x, y), 

(Hi) S and T are inverses, i.e. So T(x) = x for a — a.e. x and T o S(y) = y for v — a.e. 

y, 
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(iv) 7o = {id, T)#a is a minimiser of the relaxed optimal transport problem \3. 

(v) For J(CT,iy) (/, g) defined by h3.2S\) . the following equality holds: 

sup J(a,y){I,9)^ inf /[7]=_inf /„[T]. 

4 The main existence result in dual space 

4.1 Statement 

We give the statement of the existence result for the equations (I3.14[) - (I3.19|) in dual variables. 
This theorem is given in [SI Theorem 5.5], where a proof is obtained by an approximation 
argument. 

Theorem 4.1. Let 1 < r < oo and Uq G L^^Aq) be an initial potential density with support 
in Aq, where Aq is a bounded open set in M.^ with Aq C x [6, 1/6] for some Q < S < 1. 
Let rt be an open bounded convex set in R^. Assume that c(-, •) is given by i3.13\) and 
that (j) satisfies \2.1\) f\yi). Then the system of semi-geostrophic equations in dual variables 
\8.14-^ - f3.19\) has a stable weak solution ((7,T) such that, with v{t,-) = T{t, •)#a{t, ■) and 
w as in liS.15]} . 

(^) 

ei'-((0,r)x A), HK^, •)IIl.(a) «S II'^o(-)IIl.(a) , Vte[0,r], 
a{t, •) e W'^in), \\a{t, •) ^C = C{n, a, c(-, ■),k, K{), V t G [0, t], 

(Hi) 

l|w(t,-)llLoo(A) s;C = C(f7, A, V< e [0,t], 

where A is a bounded open domain in R'^ containing supp{v), such that A C x [(5, j-] for 
some < (5 < 1. 

4.2 Alternative proof of Theorem 14.11 

Ambrosio and Gangbo in [2J proved a general result on the solution of Hamiltonian ODEs 
that may be applied to the dual space existence problem. We briefly summarise the work in 
[2] and then describe how it can be applied to the compressible semi-geostrophic equations 
we have been considering. 

Remark 4.2. In what follows we deal with concave rather than convex functions. Hence 
we replace all definitions in [5] regarding subdifferentiability and (A)— convexity with the 
following definitions regarding superdifferentiability and (A)— concavity. This replacement 
does not affect the results of [2^. 

Definition 4.1. Let H : P^^(R^) — >■ (-co, +cxd] be a proper, upper semicontinuous function 
and let v G D{H). We say that v G L^(t^;R'^) belongs to the Frechet superdifferential dH{y) 

H(yh) < H{u) + f v(y) • (R^;'' (y) - y) iy{y) dy + o{W2 [i^, i^h)) 

as Vfi — i- V. We denote by doH{v) the element of dH{v) of minimal L'^(i/;M.^)~norm. Note 
that, by the minimality of its norm, dQH{v) belongs to dH{v) nT^P^g(R'^). 
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In the following lemma we state a continuity property of optimal plans or maps. 

Lemma 4.3. Assume that {vn}'^=i, {A*n}^i o.fe bounded sequences in P^^{M.^) narrowly 
converging to v and ^ respectively. Assume that ro(i^, /x) contains a unique plan 70 induced 
by the optimal map Rf^ : R-^ — ^ M'^ (see Remark [KJ\) . Then 

lim / giy,R^,:{y)yn{y)dy= f g{y,Il^Ay)Hy)dy, 

where R(^" is optimal in the transport of Vn to pin, o,nd for any continuous Junction g : 
X A satisfying 

lim My^j^ = 0. 

|(y,x)K + oo |y|2 + |x|2 

Assume furthermore that there exists a closed ball Br, of finite radius r, containing the 
supports of pin and fi. Then there exist Lipschitz, convex functions Un,u : M'^ — >■ MU{+oo} 
such that Vu„ — R(^j^ Vn — ci-e. in M'^ and Vm = R(^ v — a.e. in M.^ . In addition, there exists 
a subsequence {nk}'^i of integers such that 

Vw„fc -> Vu a.e. in M^. (4.1) 

Proof. See [3J Proposition 7.1.3] for proof. □ 

Definition 4.2. Let H : P^cO^^) (-00, +00] be proper and let A £ M. Let tti : x R3 : 

(x, y) — s> X and 7r2 : R'^ x R'^ : (x, y) — s> y 6e the first and second projections o/R"^ x R"^ onto 
M.^ . We say that H is X— concave if for every vi, V2 £ P^^{M.^) and every optimal transport 
plan 7 G ro(;^i, 1^2) we have 

Hiiy^t)) ^ (1 - t)Hiiyi) + tH{y2) - \t{l - t)Wi{vi, V2) 
for all t £ [0, 1], where V{^t) = ((1 ^ O""! + ^'^2)#7- 

Proposition 4.4. Let H : Pl^iM?') ( —00, +00] be upper semicontinuous and X— concave 
for some A £ R and let v £ D{H). Then, the following condition is equivalent tow G dHiv): 

H{vn) ^ H(i,) + / v(y) • (R^;'' (y) - y)z.(y) dy + ^Wi{v, vn) 

for alli^hePliR^). 

Proof. See [3J Proposition 4.2] for proof. □ 

Following [2], we define Hamiltonian ODE's as follows: 

Definition 4.3. Let H : P^^(R'^) — > (— oo,+cxd] be a proper, upper semicontinuous function. 
Define the linear transformation J : R'^ — >■ R'^ by 

J{My),v2{y),y3{y)) = f;M-My),yiiy),o), (4-2) 

for all v(y) £ R'^. We say that an absolutely continuous curve v^i^ : [0, r] — >■ D(H) is a 
Hamiltonian ODE relative to H , starting from vq £ P^^iM?), if there exist v^t) £ -^^(i^(t); R'^) 
with ||^2(j^ ) £ -^^(0, t), such that 

1 V(,) £ r,j„p2^(R3) n dH{v^t)) for a.e. t. 
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The main result of concerns Hamiltonians H satisfying the foUowing properties: 

(HI) There exist constants Cq G (0,+oo), Rq G (0, +cx)] such that for all v G Pqj,(M'^) 
with Wiiy-iVQ) < Rq we have v G D{H), dH{i') =/= and v ~ daH(y) satisfies 
|v(z)| < Co(1 + |z|) for t/-a.e. zgR3. 

(H2) // V, Vn G P^j,(R'^), sup„ W2(i^n, t'o) < ^0 md Vn ^ V narrowly, then there 
exists a subsequence n[k) and functions -Vk, v such that Vfc = dfjH{i'ri{k)) ^n{k) ~ Ci.e., 
V = doH{iy) V — a.e. and v^^ — )• v a.e. in M"^ as k — )• +cx). 

To ensure the constancy of H along the solutions of the Hamiltonian system we consider 
also: 

(H3) H : P^|,(M'^) — > (—00, +00] is proper, upper semicontinuous and X— concave for some 
A G K. 

For Hamiltonians H as above, the following result holds (see [2} Theorem 6.6] for full 
details and proof): 

Theorem 4.5. Assume that (HI) and (H2) hold for H{v) and that t > satisfies 

Cqt ^2^{l + e(25Co^+i)-(l + M2{vq))) < Ro- (4.4) 

Then there exists an absolutely continuous Hamiltonian flow Vf^i) G P'^^{E?), v^^i^ : [0,t] 
D{H) starting from G P^^{M.^), satisfying i4-^ > such that the velocity field V(^t) coincides 
with doH{i'(^t)) for a.e. t G [0, r] . Furthermore, the function t is Lipschitz continuous. 

Finally, there exists a function l{r) depending only on r and Cq such that 

Vq ^ nir a.e. on Br for all r > => i>^t) ^ o.e. on Br for all r > (4.5) 

and 

Vq ^ Mr a.e. on Br for all r > ^ Mi(^r) o.e. on Br for all r > 0. (4.6) 

// in addition (H3) holds, then t M> H{v(^f^) is constant. 

We now apply Theorem 14.51 to our problem in order to obtain directly the existence 
result of Theorem 14.11 

Let A, 57 be as in (13.201) . For v G Pac(A), define the Hamiltonian H by 

H{v) = inf [£{v, ti)+K, f {^i{^)r dx j , (4.7) 

where 

£{u,^l)^_m^ /5(y,S(y))Ky)dy (4.8) 

with £(y,x) = c(x, y) defined by p.l3p . We begin with the following: 

Proposition 4.6. Let ft and A satisfy i3.20p . Let the Hamiltonian H{v) on Pac(A) he 
defined by |^.7| j. Then H is superdifjerentiable, upper semicontinuous and {~2)-~ concave. 
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Proof. Given i> G Pac(A), denote by a the minimiser in (|4.7p . The existence and uniqueness 
of this minimiser foUows from Theorem 13.31 For any G Pqc(A) we have 

i^K) - inf + /(M(x))"dxl (4.9) 



Jn 



First, recall that we can guarantee that there exists a unique optimal transport map Rj^'' 
from u to with respect to the Wasserstein cost function d{y,yh) = — y/ip; see 
Remark O 

We consider now transport with respect to the cost function c(y, x) = c(x, y) given by 
p.l3p . Let be the optimal map in the transport of i/ to a and let S^^ be the optimal 
map in the transport of Vh to a. Therefore, we have 

inf / £(y, S(y))Ky) dy ^ f £(y, S-(y))Ky) dy 

and 

inf / £{y,S{y))My)dy = f c{y,S:jy))My)dy. 

The existence of and S^^ follows from TheoremlXSl Note that, since (S^o(R|^'' = 
a and since S^^^ is optimal in the transport of Vh to cr, we have 

/ c{y,S:^{y)),.^{y)dy^ f 5(y, SXRj;'-)-i(y));.„(y) dy. 

J A J A 

It follows that 

- H{iy) < Sii^h, a) + Ki /(c7(x))'' dx ~ a) - Ki /(a(x))'' dx 
Jn Jn 

= [ c{y,s:jy))My)dy- f a(y, (y))Ky) dy 

J A J A 

< / c~(y,SXR;^'^)-i(y));.„(y)dy- / 5(y, S-(y))i.(y) dy 

J A J A 

c(Rj:'^(y),S^(y))Ky)dy- / £(y,S^(y))Ky)dy 

V J A 

din:- (y), S^(y)) - c~(y, (y))] u{y) dy, 



[ Vc(y, S^(y)) ■ [R;^-(y) - yHy) dy + o{W2{i^, vn)). 

J A 



(4.10) 

where we have used (|2.5p . 

Hence, using Definition 14. 1[ we conclude that Vc(y, S^(y)) S dH{v). Thus, dH{v) is 
non-empty, H is superdifferentiable and we can use [14, Proposition 10.12] to conclude that 
H is semi-concave, i.e. 

H is (-2) - concave. (4.11) 

Also, from the narrow continuity of f (•, •) (see |S1 Theorem 3.4]) and the uniform conver- 
gence of a as the minimiser of (|4.7p (see jHl Lemma 4.3]), we have that 

H is upper semicontinuous. (4.12) 
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From ((ilT|) and we have that (H3) holds. □ 

The foUowing proposition yields an alternative proof of the main result of [8] ; Theorem 

Proposition 4.7. Let 1 < r < oo and vq G L''(Ao) be an initial potential density with 
support in Aq, where Aq is a bounded open set in with Aq C x 1/(5] for some 
< (5 < 1. Let n and A satisfy i3.20\) and let the Hamiltonian H be defined by ^-T^ - 
Then, there exists a Hamiltonian flow i>(^t) G ^ac(A) and constant r > such that 

^l^(t) + V ■ (J(v(4))i^(4)) = 0, 1/(0) = i/Q, te{0,T) 

where J(v(-j')) = w a.e. in [0, r]. 

Proof. We compute dQH{v) (as defined in Definition UT]) explicitly to show that the condi- 
tions required to apply Theorem 14.51 hold. From the definition of J in (14.21) , velocity fields 
transporting v will have vanishing components in the 7/3 direction so that we need only 
consider variations of v in the (j/i, y2)~directions. Thus, to characterise the elements of 
dH{v), we let Lp e Cj?°(R'^) be arbitrary, but with no j/3 dependence, and set 

gs(y) = ((5.)i(y), (g.)2(y), (5.)3(y)) = y + sV^{y). 

Note that {gs)d.{y) = 2/3 and, for \s\ sufficiently small, gs is the gradient of a convex function, 
since gs(y) = V(^y^ + s(f). Define Vs = Ssf^>^- Denote by as the minimiser in 



H{iys)= inf U{iys,t^)+Ki {fi{^))^dii 

The existence and uniqueness of the minimiser CTs follows from Theorem lB.BI Let ^ G dH{v). 
Combining the (—2)— concavity of H and (I4.12p with Proposition 14. 4[ we obtain 

H{us) - H{u) - f e(y) • (R^^ (y) - yHy) dy + Wi[y, y^) ^ 0. (4.13) 

J A 

Brenier's polar factorisation theorem states that any suitable mapping from v to can be 
uniquely factorised as the composition of a measure-preserving mapping and the gradient 
of a convex function (sec, for example, [13, Chapter 3]). Since, for |s| sufficiently small, gs 
is the gradient of a convex function, we conclude that 

Wi{v,vs)= I \y-'K'{y)?^{y)dy= I \y - Qs{y)?y{y)dy ^ [ \V^iy)\^Hy)dy 

JA J A J A 

and 

/ e(y) • (R^^ (y) - y)^(y) dy = ( ^y) ■ (gs(y) - y)^^(y) dy = s [ ^y) ■ ^^{yMy) dy. 

JA JA JA 

Combining this with (j4.13p . we therefore obtain 

/ ^(y) • V'piyHy) dy + f |V(p(y)|2i.(y) dy ^ H{y) - H{vs) 

J A J A 

^ £{iy, as) - £{iys,o-s) 

= / S(y,S^iy))i/(y)dy- / 5(y, S^^ (y))z.,(y) dy 

J A J A 

^ / c(y, S^; o g,(y))z.(y) dy - [ 5(y, S^; (y))i^.(y) dy 

JA JA 
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/ £(g7i(y),S-=(y))i.,(y)dy- / 5(y, S^; (y))i.. (y) dy, 

J A J A 



(4.14) 



since gs^i^ = Vg. Here SJ^'' denotes the optimal transport map from v to and SJ^= 
denotes the optimal transport map from Va to with respect to the cost function c(-, ■). 
The existence of SJ^=^ and follows from Theorem 13.51 
Note that 

gs '(y) = y - sV(/j(y) + y VV(y)V(^(y) + e(s,y), 

where e is a function such that |e(s,y)| ^ |sp||<y9||c'3(R3). 

Combining this expression for g^"^ with (|4.14p and using -^"-P = 0, we conclude that 

-s /e(y) • V¥.(y)Ky) dy + / |V(p(y)PKy)rfy / [5(g7i(y),S^;(y)) -5(y,S:;(y))] i/,(y)dy 

J A J A 

^/L{|(.9.)r'(y) - {SZi\{y)\^ + \{9s)-^\y) (5-;)2(y)r} + 0(S^;(y)) 



2/3 



^/clr{|yi - (5.^;)i(y)r + 1^2 - (^.^;)2(y)r} + 0(s^;:(y)) 



1 

A 2/3 



12 cor 



y3 



^/L {I2/1 - (5.^;)i(y)|' + 12/2 - (^:;)2(y)r} 
s^:(y)-y 



i^s{y) dy + o{s) 



2/3 



V(p(y)j^s(y) dy + o{s). 



Recall now that ^ in ^ac(A) and CTs ^ cr in Pac(^^) as s — > 0, hence Theorem 14.31 gives 

's^(y)-y^ 



s / ^(y)-V(^(y)Ky)c?y+s' / |V(^(y)|^i^(y)dy s/, 



/ cor 



2/3 



•V'^(y)j^(y) dy+o{s). 



Dividing both sides first by s > 0, then by s < and letting |s| — > we obtain 

's^(y)-y^ 



e(y)-v^(y)Ky)c?y = /< 



2/3 



V(y5(y) i^(y) dy. 



Thus, we have that 



Ji^.ay)) = •^(/cor(^^^)), where tt. : L2(^.a) ^ T.p2^(A) 
denotes the canonical orthogonal projection, where the tangent space is defined in 
The minimality of the norm of doH then gives 



c2 

' cor 



s^(y)-y 



2/3 



w(y), 



(4.15) 



where w is defined as in p.l5|) . 

We can now check directly that conditions (HI) and (H2) hold. Using Theorem 13. 4l and 
Theorem 13.51 we may conclude from (I4.15P that 



J(aoi^(;^)) =w(y) = J(V5o(y)). 
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Condition (HI) then follows from the fact that, as a solution of the dual problem p.23p . 
go G W^'°°{A). Condition (H2) follows from the stability of optimal maps (see, for example, 
[HI Corollary 5.23]). Hence we may apply the result of Theorem 14.51 to conclude that there 
exists a Hamiltonian flow such that 

^i^(t) + V ■ (J(v(t))j/(t)) = 0, 1/(0) = i/Q, ie(0,T) 

where J{'V(t)) — ^(<9o-ff (^^(t))) for a.e. t e [0, r]. By (|4.15p . this then completes the proof 
that the dual space continuity equation (|3.14p . with velocity field defined as in (I3.15p . is 
satisfied. In addition, from (H3) and the definition of J, the energy associated with the 
flow is conserved. 

□ 

Concluding remarks on the dual space result 

The result of the previous section yields a proof of Theorem 14.11 alternative to the existing 
proof of [8 . Indeed, by the deflnition of w in Proposition 14.71 we have that (ct, T) is a 
stable solution of (P?Til) - ((XTOl) . where T = (see TheoreminU]). Theorem |3TT](i) follows 
from (j4.5p . (14.61) : Theorem 14. II (ii) follows from the fact that a is a minimiser of (j4.7D (see 
[8j Theorem 4.2]); Theorem 14.11 (Hi) follows from the deflnition of w in terms of Vgo and 
the fact that, as a solution of the dual problem p.23p . go £ W^'°°{A). 

5 Lagrangian statement of the equations in physical 
space 

We aim to prove the existence of a weak Lagrangian solution of the fully compressible 
semi-geostrophic system (I3.1l) - (|3.6p . We will show that, with additional regularity, a weak 
Lagrangian solution would determine a weak Eulerian solution. 

We begin by rewriting the system (I3.ip - p.6p in a form that enables us to state the 
equations in Lagrangian form. These new equations (with prescribed initial conditions) are 

AT(t, x) = /cores X [T(t, x) - x] , (5.1) 

dta{t, x) + V • {(T{t, x)u(t, x)) = 0, (5.2) 

Vxc(x,T(t,x)) + KifiV((a(i,x))''-i) = 0, (5.3) 

u-n==0 on[0,T)xdn, (5.4) 

cr(0,x) = (To(x), T(0,x)=To, To #(70 e i''(Ao), Ao C compact, (5.5) 

where r G (1, oo), c(x, T(t, x)) is deflned in (j3.13p and Ki is defined as in (j3.9p . 
Proposition 5.1. A solution of l5.1\} - ![575\} determines a solution of the original system 
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Proof. Given a solution (u,T, cr) of (|5.ip - (l5.5p . we set the function u in (|3.ip - p.6p to be 
equal to the function u in (|5.ip - (|5.5l) and we define 

0(i,x) :=r3(i,x), p(i,x):=5j^, u5(i,x):=/eore3X [T(i,x)-x]. (5.6) 

Then p.2p follows from (|5.6|) and equation DtT^ — of (|5.ip . Using this together with 
(|5.6D and the fact that (j5.2p is satisfied, we see that p.3p holds. 

In order to show that (|3.4p holds, we first rearrange p.5p to obtain p = R'^pl~f'^{p9)'^ . 
Thus, 

Vp = i?''pi^l''«(p0)'^~iV(p0). (5.7) 
Then, substituting (|5.6I) into (j5.3p gives 

/cores X uf + V0 + 6'Kii:iV(p6i)''"i = 0. 



Recalling that Ki — Cy ( — ^ J , this becomes 

\ ^rcf / 



/cores X u9 + V(/) + enc, ( —] VipO)^-^ = 0, 

VPrcf/ 

and therefore 

/cores X + V0 + -KiK - l)c, f — ) (p0)''-i V(p0) = 0. 

P VPrcf/ 

Then, using ^^(xiv), (^Al (xv) and ((5J)) we obtain 

/cores X uf + V(/) + = 
P 

and thus conclude that p.4p is satisfied. 

Finally, we obtain p.ip using p.4p together with the first two components of (j5.ip 
and the definition of u^' in (15.61) . Hence, we obtain a solution of the original system of 
semi-geostrophic equations. □ 

We define a weak solution of (|5.ip - (|5.5p as follows: 

Definition 5.1. Let fl be an open bounded convex set in M^. Let u : [0, r) x — > R'^ 

satisfy u e L'^{[0,t) x n) and let T : [0, r) x rj -J> satisfy T e L°°([0,r) x fi). ^/so, 
fei a{t,-) e VF^'°°(f2) /or a// t e [0,t]. T/ick (u,T, ct) is a weak (Eulerian) solution of 
(EjP-fSfaj 2/ 

/ {T(i, x) • [dMt, x) + (u(i, x) • V)(/)(t, x)] + /co,es x [T(i, x) - x] • ^{t, x)}a{t, x) dtdx 

J[0,r)xn 

+ / To(x) • (^(0, x)ao(x) dx = 0, (5.8) 
Jn 

for any ip G Cc([0,t) x 17), and 

{9f(/'(t, x) + {u{t, x) • V)-^^*, x)}CT(i, x) d<dx + / V(0, x)cro(x) dx = 0, (5.9) 



'[0,t)xJ2 

/or any V e Ci([0,r) x 17). 
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Given a solution T{t, x) and x(t, y) = S{t, y) of the system (|3.14p - (|3.19p in dual coordi- 
nates, formally we have that (u, T, a) satisfy (|5.ip - (|5.5p . However, due to the low regularity 
of T{t, x) and thus S{t, y) obtained as a weak solution of p.l4p - p.l9p . this argument is not 
rigorous. Indeed, the regularity obtained in [8] may not be enough to provide a solution of 
the Eulerian problem. 

Remark 5.2. In a recent development, [?], the existence of weak Eulerian solutions of 
the incompressible semi-geostrophic equations on the 2-dimensional torus has been proved. 
Whilst the authors expect that this result could be extended to three dimensions, the 
question of whether it is possible to prove that weak Eulerian solutions exist for the fully 
compressible semi-geostrophic equations remains open. 

Instead, we seek to find weak Lagrangian solutions of the problem (j5.ip - (l5.5l) . Here we 
define such solutions and then we use the methods of 6, to prove their existence. We define 
the Lagrangian flow map 

F :[o,t) xn^n 

corresponding to the velocity u. Note that F maps 57 to itself so that the boundary 
conditions of the problem are respected. Then, we can rewrite the system (j5.ip - (l5.5p in 
terms of this Lagrangian flow F and define the corresponding weak solution (F, T, cr). 

Definition 5.2. Let 57 be an open bounded convex set in and let t > 0. Let T € 

L°°{[0,t) X n) and let a{t,-) S W^'°°{n) for all t € [0,t]. For I q < oo, let F : 
[0, r) X — > 51 &e a Borel map, satisfying 

F e C([0,r);L9(f];M3)). (5.10) 

Then (F,T, cr) is called a weak Lagrangian solution of lid. l]) - ^i5~^) in [0, r) x 57 if: 

(i) For (To ~ x e 57, 

F(0,x)=x, T(0,x) =To(x). 

(ii) For any t > Q, the mapping 

F(t) =F(t,-) : [0,t) X 51 ^57 
pushes forward the probability measure cro(') to a{t, ■), i.e. 

F(t)#ao(-)=^(t,-). (5.11) 

(Hi) There exists a Borel map 

F* : [0, r) X 57 ^ 57 

such that, for every t S (0, r), the map 

F^t) = F*(i,-) : 57 51 

satisfies 

F(t) ° F(t) W = ^ ^"'^ ° W = X' 
for (7 — a.e. x G 57, and 

F*#a{t,-)=aoi-). 
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(iv) The function 

Z(t,x) :==T(t,F(,)(x)) (5.12) 

is a weak solution of 



dtZ{t, x) ^ f,ore3 X [Z(t, x) - F(t, x)] m [0, r) x 17, 
Z(0,x)=To(x) inn, 



(5.13) 



m i/ie following sense: 
for any ip G C^{[0,t) x f2; 



/" {Z(t, x) • dMt, x) + fcores X [Z(t, x) - F{t, x)] • ^(t, x)}ao(x) dtdx 

J[o,T}xn 

+ [ [To(x) •¥.(0,x)]f7o(x)dx = 0. 



(5.14) 



Remark 5.3. Let us comment on Definition 

• Continuity in time of F considered as a map on [0, r) with values in i*(f2) as required 
in (jS.lOL combined with the initial condition for F in (i), implies that 

||F(,)-m^,(^^=0. (5.15) 

Furthermore, the continuity property (I5.10p may be interpreted as continuity of par- 
ticle paths in physical space. 

• Property (ii) is the Lagrangian form of the mass conservation principle in equation 
(|5.2p with boundary condition ([57 



• Equation (|5.13l) is the Lagrangian form of (|5.ip . 

• Equation (|5.14p is derived by multiplying (I5.13P by ao{x) as well as the test function 
and then integrating by parts, as in [IT, Definition 3.4]. The reason for this choice will 
become clear in the proof of Proposition 15.41 

• We have omitted equation (15. 3p since this holds as a result of the energy minimisation; 
see [H Theorem 4.2]. 



In order to justify Definition 15.21 we must show that a weak Lagrangian solution corre- 
sponds to a weak (Eulerian) solution of (j5.ip - (j5.5p . as defined by Definition [5Tj Indeed, we 
prove that, with additional regularity property dtF G L°°([0,t) x fl), a weak Lagrangian 
solution (F, T, a) as defined in Definition 15 . 2 1 determines a weak (Eulerian) solution of (|5.ip - 
(15.51) and, furthermore, that a smooth Lagrangian solution determines a classical solution 
of (j5.ip - (j5.5p . This is the content of the following result: 

Proposition 5.4. Let be an open bounded convex set in and let r > 0. Let (F, T, a) 
be a weak Lagrangian solution of i5.1\) - ![57S\) in [0,t) x fl. 

(i) If dtF e L°°([0,r) x n;R^), then the function 

u(i,x) :-(aiF)(t,F^,)(x)) (5.16) 

satisfies u G L°°{[0, r) x il; M."^) and (u, T, a) is a weak Eulerian solution of i5.1\) - ![575\) 
in [0,t) X n in the sense of Definition \5.i\ 
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(a) //(r,F*,T) e C'^{[0,T]xn), then the function fO^j satisfies u <E C^{[0,T]xn;M.^), 
and (u, T,ct) is a classical solution of L5.1\) - l[575\) in [0,t) x fl. 

Proof. Let us first prove (i). Since F* is a Borel map and, by our additional regularity 
assumption, dtF G L°°([0, r) x 51), we know that the right-hand side of (|5.16p is a bounded 
measurable function. Therefore, u e L°°{[0,t) x ft). Now, in order to prove that (u,T, cr) 
is a weak Eulerian solution of (j5.1l) - (j5.5p . we must show that (I5.8P and (j5.9p hold. We begin 
with (|5.9p . Let -0 G Cc([Oi ''') ^ that the support of ■4> in i is a closed subset of [0,t), 

and fix t G (0, r). Note that, since a(t, •) G VKi'°°(fl) for all t G [0, r] (by TheoremO (li)) 
we have that (T{t, •) G i^([0,r) x fJ) and, since F(t)#(To(-) = •), we can apply (|2.5p to 
yield 

(at7/.)(i,F(t)(x))ao(x)dx= / atV^(i,X)a(t,X)dX. 
:! Jn 

Then, applying the chain rule, integrating both sides with respect to t and using our 
assumption that dtF G L°°([0,r) x fi), we obtain 

{9t[V'(t,F(t)(x))] -5tF(4)(x) . (V0)(i,F(4)(x))}ao(x)dtdx 

[0,r)xO 

= [ dt^{t,X)a{t,X) dtdX. 
Now, using (I5.15P and the fact that iP{t,-) = (by compact support), we get 
9t[V'(i,F(t)(x))](To(x)dtdx - / ip{0,x.)ao{x.) dx. 



[o,T)xn 



Hence, 



V'(0,x)(7o(x)dx- / atF(t)(x)-(VV')(t,F(t)(x))(7o(x)dtdx 

[0,T)yA1 

dt^{t,X)a{t,X)dtdX. 

[0,t)xO 

If we make the change of variables X = F(j)(x) in the second integral above, then, by (ii), 
(Hi) in Definition 15.21 and by (|2.5|) . we have that 

- / V(0,x)ao(x)dx- / {dtF(t)){Fl,.{X))-Vi^{t,X)a{t,X)dtdX 

Jn J[0,r)xO 

dtip{t,X)a{t,X)dtdX. 

'[0,r)xO 

Then, rearranging and using the definition of u in (j5.16p . we obtain 

/ {dt^it, X) + u{t, X) • \7^{t, X)}a{t, X) dtdX + [ V(0, x)cro(x) dx ^ 0. 

J[o,T)xn Jn 

Changing notations X to x gives us (|5.9p . 

We now prove that (|5.8I) also holds. By the properties of F and T in Definition 15.21 
we have that Z(i,x) as defined in (I5.12p satisfies Z G L°°{[0,t] x fl). Also, applying the 
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definition of Z(t,x) in (|5.12p to equation (|5.14p gives 

/ {T(t,F(t)(x))-at^(t,x) + /cores X [T(t,F(t)(x)) -F(t)(x)] •(p(i,x)}ao(x)dtdx 
+ [ [To(x) •(p(0,x)](7o(x)dx = 0, 

(5.17) 

for any ip e ([0, r) x fi). Now, since is a bounded set and F(4)#cro(') — /^{t, •) for all 
t S [0,t), equation (|2.5p allows us to make the change of variables X = F(4)(x) in the first 
integral of (|5.17l) . Thus, by (in) of Definition 15.21 we have that x = F*j^(X) for a — a.e. 
X G for every t G [0, r) and then, from (|5.17p . we obtain 

/ {T{t, X) • dMt, (X)) + /cores X [T{t, X) - X] • p{t, FT. {X))}a{t, X) dtdX 
j[o,T)xn 

+ / [To(x)-¥.(0,x)]f7o(x)dx = 0, 

(5.18) 

for any (p G ([0, r) x 57). We now show that (15.181) also holds for all (p such that 

¥.eL°°([0,r) xf]), 
9t^ei°°([0,T) X 17), (5.19) 
supp{if) c[0, T — e] X 17 for some e > 

In order to do this, we construct an approximating sequence for such ip. Let us extend <p to 
(— oo, oo) X 17 by defining, for x e 17, ip[t, x) = <p(— x) for i < and 'p>{t, •) = for t ^ r. 
Then, we let ft. > and define 17/i = {x G 17 : dist{x.,dil) > h}, where 917 denotes the 
boundary of H.. Thus, tpx^^ ^'^'^ defined on x R^, where xn^ denotes the characteristic 
function of the set ^h- Next, let jh{t,x) = 7^i(^%f^)j where j(-) is a standard mollifier, 
and let A: > i be an integer. We then have that functions ipk = {'PXo,4,^)*jh, with h = ^ < e, 
satisfy 

p>k e CI{[0,t) X n) with \\ipk,dt(pk\\L'=--uo,^}xn) ^ C, 
where C does not depend on fc, and 

{ipk,dtPk) ^ iVidt^) a.e. on [0,t) x 17 as fc — ?> oo. 

Thus, by the dominated convergence theorem, 

lim / {ipk,dtipk){t,:>i)(7o{^)dtdx= / {ip, dt(p) {t, x)<to {x) dtdx. 
Also, since ip e L°°([0,t) x n) by ((5l9| and F*^^#a{t, •) = cto(-) , we have from ([231) that 



((pfc,5f.^fc)(i>x)CTo(x)dtdx = / {lpk,^t^Pk)it,F*^•^(X))(7{t,X.)dtd'X 

[0,r)xO J[0,T)xa 



and 

/ 



((^,9t^)(i,x)ao(x)dtdx= / (^,5t^)(i,F^,)(X))a(t,X)dtdX. 

:"2 Jfo.r'ixn 
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Hence, 

lini / {^k,dt^k){t,F%JX))a{t,X)dtdX= f F^^ (X))a(t, X) dtdX. 

J[o,T)xn J[o,T)-xn 

Then, since ft is bounded, T G L°°{[0, r) x ft) and (j5.18|) holds for each (pk, it follows that 
([ET5|) holds for ip satisfying ((5l^ . Now, let 

^(t,x)=77(t,F(t)(x)), (5.20) 

where 77 G C;^([0,t) x il). Then, since supp{ri) d [0,t ~ e] x ft and 9tF e L°°([0,t) x fi) 
by our additional regularity property, we have that (p satisfies the conditions in (|5.19p . 
Therefore, (|5.18p holds for ip as defined in (|5.20p . We also see from (I5.20p and use of the 
chain rule that 

dMt,^) = (9*r;)(i,F(,)(x)) +atF(,)(x) • (Vr;)(i,F(,)(x)). 
Thus, by property {in) of Definition 15.21 we have 

(^{t,F^,)(X))=,7(t,X) 

and 

dMt,Fl,^m) = dMt,X) + [(9tF)(i,F^,)(X))] . Vvit,X) 

for cr — a.e. X G for every t G [0,t). Then, inserting p{t,x.) = 77(t, F(t)(x)) into (|5.18p 
gives 

/ {T{t, X)-[dtv{t, X)+[{dtF){t, F%. {X))]-Vr,{t, X)]+/eo.e3 x [T(t, X)-X].r?(t, X)}a{t, X) dtdX 
J[o,T)xn 

+ [ To(x)-77(0,F(o)(x))ao(x)dx = 0, 

and, using the definition of u in (|5.16p together with property (i) of Definition 15. 2[ we 
obtain 

/ {Tit, X)-[dtvit, X) + {u{t, X). V)77(i, X)]+/eore3 X [T(i, X)-X]-7j{t, X)}a{t, X) dtdX 

J[0,t)xS1 

+ / To(x) •77(0,x)ao(x)dx = 0. 

Finally, changing notations X to x and 77 to ip gives (|5.8p . proving statement (z). Then, 
statement (ii) follows directly from (i). □ 

We now state the main result, which we will prove in Sections [6] and [7j 

Theorem 5.5. LetQcM.^ be an open bounded convex set. Assume that 

To#<7o ei''(Ao) 

for r G {1,00), where Aq C R"^ is compact. Then, for any t > 0, there exists a weak 
Lagrangian solution (F,T, cr) of S5. ip -i TO)) in [0,t) x Q, where T G L°°{[0,t) X fl) and 
Ii5.10\} is satisfied for any q G [l,cx)). Moreover, the function Z(t,x) defined by 15.12\) 
satisfies Z(-,x) G VF^'°°([0, r)) for (Tq — a.e. x G fJ, and i5.1S\) is satisfied, in addition to 
the weak form ^5.14^ , in the following sense: 

dt'Zi{t,:x.) = fcor^s X [Z(t,x) — F(t,x)] for uo — a.e. x G 51 and every t G [0,t), 
Z(0, x) = To(x) for (To - a.e. x e fl. 
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To prove this theorem, we show the existence of a Lagrangian flow map $ in dual 
variables. We then define F in terms of $. 

Throughout the next two sections, we will use tr and v to denote the measures cr(i, •) 
and v(t^ •) considered at some fixed time t £ [0, r). Thus, if we write v — a.e. y G M.'^, for 
example, then we mean v(t, •) — a.e. y G M"^ for a fixed t G [0, r). 

6 Existence of the Lagrangian flow map in dual space 

In this section we use the results of [T to show existence of the Lagrangian fiow map in 
dual space. 

Let f2, A be as in Theorem 14. 11 let Tq, (Jq be as in (|5.5|) and let vq = To#cro- Then, by 
Theorem 14. II there exists a solution (T, S,(t) of the system p.l4p - p.l9p . with initial data 
vq, satisfying all assertions of Theorem 14.11 We use Theorem 13. 4[ Theorem 13.51 and the 
definition of w in (j3.15p to see that the dual space velocity can be written as 

w = /coJysea X Vgoly) (6.1) 

where go is a solution of the Kantorovich dual problem (j3.23l) . Thus, the vector field w is 
divergence free. Furthermore, since the potential temperature 2/3 is assumed to be bounded, 
we can use (|6.ip to obtain the following properties of w: 

• by the semi-concavity of go in A, it follows that w(i, •) G BViod'^^) f^i' '^■S- t € (0, t); 

• by Theorem Q] we have that w G L°°{[Q,t) x A). 

Since u has compact support in [0, r] x M"^, we can modify w away from A with a standard 
smooth cut off function C G C°°(R'^) so that the modified function w := satisfies 

w G L°°([0,r) X M^)^ w(i, •) G BV;oc(M^) for a.e. i G (0,r), (6.2) 

V • w(t, ■) = Q in for every t G [0, r). (6.3) 

and 

w = w in A. (6-4) 

These conditions enable us to apply the theory of JLj to the transport equation p.l4p with 
w replaced by our modified velocity w: 

Lemma 6.1. There exists a unique locally bounded Borel measurable map $ : [0, r) x K"^ — 
M.^ satisfying 

(i) $(-,y) eW^-°°{[Q,T)) for vo- a.e. y G M^- 

(ii) $(0, y) = y for vq - a.e. y G K^; 

(Hi) for — a.e. y G M'^, 

9t<i>(t,y)=w(t,$(i,y)); (6.5) 

(iv) there exists a Borel map $* : [0,t) x R-^ — > R-^ such that, for every t G (0,r), the 
map <i>^jj : R"^ — > M'^ is Lehesgue-measure preserving, and such that <^*^^-^ o (y) = 

*(t) ° *(t)(y) =y forv- a.e. y G R^■ 

(v) <&(<, •) : M'' — > M.^ is a Lebesgue-measure preserving map for every t G [0, r). 
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Proof. The proof is essentially identical to that of [51 Lemma 2.8]. 

□ 

We now show that the image of the flow map $ is contained in A, and therefore corre- 
sponds to the velocity field w. 

Lemma 6.2. Let A be as in Theorem \4.1\ Let $ be the map defined in Lemma \6.1\ and let 
w be defined as in 113.15]) . Then 

^{t,y) C A forvQ-a.e. y e Toi^l) and every t e [0,t). (6.6) 

In particular, 

aj$(i,y) ^ w(t, $(i,y)) forPQ-a.e. y e To{n) and every t e [0,t). (6.7) 

Proof. Firstly note that, since To^cto — '^o and the measures ctq and z/q are contained 
within n and A respectively, we may assume 

To : 17 -> A. (6.8) 

Therefore, for (Tq ~ ^-e. x G f7, we have that To(x) G A. 

We begin with the vertical component of $3. From (|6.4I) and (|6.5I) . we have that 
dt^3{t,y) = W3{t,^{t,y)) = for all y G and every t G [0, r). Therefore, we have 
^sit^y) = ^3(0jy) — 2/3 for z^o— a.e. y G M"^ and every t G [0, r), where we have used 
Lemma r6. II (ii). We therefore conclude that S < $3(t,y) < | for j^o— a.e. y G To(ri) and 
every t G [0, r), where 6 is as in Remark \3A\ 

The proof that the horizontal components of $ stay inside A proceeds in a similar way 
to that of [5J Lemma 2.9]. 

□ 

Now, we wish to prove that when {i^,T) is a weak solution of p.l4p - p.l9p then v is 
a weak Lagrangian solution of the transport equation (j3.14[) . i.e. v satisfies the property 

f = $#t'o. 

First, we state the following result from [3]: 

Proposition 6.3. Let v(t, ■) be narrowly continuous Borel probability measures solving the 
continuity equation 

dtv + W ■ (vv) = 
with respect to a vector field v satisfying 

\v\i^dtdy < +00, (6.9) 

[0,r)xE3 



sup |v| + Lip{v, B) 

B 



dt < +CX3, (6.10) 



for every compact set B C M''. Then, for i^q — a.e. y G M'' the characteristic system 

*(0,y) = y, ^<i>(t,y) -v(t,$(t,y)) (6.11) 

admits a globally defined solution $(t,y) in [0,t) and 

i^{t,-)^^t,-)#M-) for allte[0,T). (6.12) 
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Proof. See [3J Proposition 8.1.8]. □ 

Note that w is not Lipschitz continuous and so does not satisfy the regularity assumption 
(j6.10p . As such, we approximate w in order to prove the foUowing 

Proposition 6.4. Let il, A, r, Tq be as in Theorem \4-.l\ and let {v, T) be the weak solution 
of !li3.14-^ - [3.19\} as constructed in Theorem \4-.l\ Let <i> be the regular Lagrangian flow of vf 
defined in Lemma \6.1\ Then, for every t G [0,t], 

v{t,-)^^{t,-)#M-)- (6.13) 

Moreover, for every t G [0, t\, 

Kt,y) = J^o(**(i,y)) forv-a.e. y eM.\ (6.14) 
where the map $* is defined in Lemma \6.1\ (yv). 

Proof. Notice first that, by (|6.2p and the compact support of w , we have 

/ \<ff\Pvdtdy <+oo, (6.15) 

where p ^ 1. 

Let (je) C C°°(M'^) be a family of strictly positive moUifiers in the y variable and set 
:= ly^je, kg = (wz/) * jg, We := — . Then, since the transport equation p.l4p is satisfied 
with w replaced by our modified velocity w, and since (j6.15D holds, we may use 3, Lemma 
8.1.9] to conclude that z/g is a continuous solution of 

dtl^e + V • (WeJ/j) — 

and that satisfy the regularity assumptions (|6.10p and the uniform integability bounds 
for all t G (0, r). Moreover, 

i>e i^, kg 'wv (6.16) 



narrowly and 



lim ||We||iP(j.^;R3) = ||w||ip(^.K3), 

£— >0+ 



for all t G (0,r). 

Thus, from Proposition 16.31 we obtain 

:/,(i,-) = $e(i,-)#f^e(0,-) 

for all t G [0,t], where is a globally defined solution of 

*(0,y)=y, ^$(i,y) =w,(t,$(i,y)). 

Thus, by Theorem 14.11 (i) and the definition of v^, we may apply (12. 5p to obtain for any, 
t > and any ip G CdR^), 

I ^{<^,{t,y))v,(0,y)dy^ [ (p(Y)z.,(t, Y) dY. 
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Passing to the limit e ^ and combining (|6.16p with T, Theorem 6.6] and the dominated 
convergence theorem, we conclude that 

/ ^mt,y))My)dy^ [ v{Y)^{t,Y)dY, (6.17) 

where $ denotes the regular Lagrangian flow relative to w (and therefore, by Lemma [6. 2 1 
relative to w in A) given by Lemma \6A\ By (j2.5l) . this implies (j6.13p . Since $ is a measure 
preserving map, we use Lemma lOl fw) to conclude that the left-hand side of (|6.17p is equal 
to 

/ ip{Y)M'^*{t,Y))dY, 

JR3 

and now (j6.17D implies (j6.14p . 

□ 



7 Lagrangian flow in physical space 

Throughout this section we will assume that 17, A, r, Tq, T, v, w, w, $ are as in Proposition 
16.41 Note that, by Theorem 14. 11 we can apply (j2.5p to ctq, a, uq, v throughout this section. 

We now perform the last step of the analysis and prove the existence of a Lagrangian 
flow F : [0, r) X — )■ in the physical space. Indeed, we define F(t) : f2 ^ for t G [0, r) 
as 

F(t) :=S(t)0$(,)oTo, (7.1) 

where Tq is as in (j5.5L S(t) is the inverse of T(f) (see Theorem 13. 5p and is the La- 
grangian flow in dual space constructed in Lemma l6. II To justify this definition, we prove 
the following lemma: 

Lemma 7.1. For any t G [0, r), the right hand side of |7.-?[ ) is defined ctq ~ a-e- The 
map F : [0, r) x — > f2 defined by ( [7._?[ ) is Borel. 

Proof. Since Tq exists and is unique ctq ~ ^-e. in $7, we have that Tq exists and is unique 
on f2 \ Nq where Nq is a Borel subset of VI with CTo[iVQ] — 0. Also, since S exists and is 
unique v — a.e. in A for every t G [0, r), we have that S exists and is unique on A \ N'^ for 
every t € [0, r), where N'^ is a Borel subset of A with j/[iV^] = 0. Then, the right-hand side 
of (TrTj) is defined for all 

X en\{N^uM), 

where 

M={Xe{n\N^) : $(t)(To(X)) GiV^}. 

Note that, from its definition, M is a Borel set. 

It remains to prove that cro[M] = for every t G [0,t) 

Fix t G [0, r). Then, using that To#o-o = and thus To#o'o = i^o for aU x G \ iVp , 
and using (j6.13p as well as Lemma l6.ll (iv). we can apply (|2.5p and compute 

ao [M^t)] = do [{x G r! \ : To(x) G %]{N^)}] (7.2) 

= / ao(x)dx= / ^o(y)rfy (7.3) 

i^iy) dy = 0. 

N2 
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Thus, we can define F : [0, r) x — ^ by (|7.ip . Then, by Lemma ISTTl F is a Borel 
mapping. □ 

It remains to prove that, if F is defined by (|7.ip . then (F,T, cr) is a weak Lagrangian 
solution of (I5.ip - (|5.5p in the sense of Definition 15.21 We begin by showing that the initial 
condition for the flow in Definition 15.21 (i) is satisfied. 

Proposition 7.2. Let F be defined as in ( |7. Then, F(0,x) = x for ctq ^ ^-e. x G fl. 

Proof. By (17. ip we have that F(o)(x) ~ S(o) o $(q) o To(x) for all x e \ iVo where Nq is a 
Borel set with cro[iVo] = 0. 

By Lemma r3. 51 there exist Borel sets TVi C 11, A^2 C A with cro[A^i] = fo[-/V2] = such 
that To exists and is unique in i7 \ Ni and S(o) exists and is unique in A \ A^2- Moreover, 
if x e f7 \ [iVi U Tg ^(7V2)], then S(o) o To(x) = x. Also, by Lemma lOfn). we have that 
<i>(o)(y) = y in A \ N^, where N3 is a Borel set with i^o[A^3] — 0- 

Therefore, F(o)(x) = x for aU x e 11 \ [Nq U iVi U Tg ^iVj U N3)]. We must now show 
that (To n Tq i(iV2 U N3)] = 0. 

We have that i^o[A^2 U N^] = 0. Then, since To#cro = 1^0, we obtain 

(To [n n To\N2 u TVs)] = / fTo(x) dx = / i/o(y) rfy = i^o[^2 u n^] = o. 

JTg ^(AfaUWa) J N2UN3 

□ 

Next, we prove that property (ii) of Definition 15.21 is satisfied. 
Proposition 7.3. For every t > 0, the map F(j-) : ^ CI as defined in J?. Jp satisfies 

F(t)#(To = (T. 

Proof. In order to prove that F(t)^(To = (t, wc must show that, for any tp G 

/ (p(F(t)(x))(To(x)dx= / (^(X)(t(X) dX. 



Then the result will follow from (12.51) . 

Let if G C(R^). From the definition of F(t) in (|7.1I) we have 

/ <^(F(t)(x))cro(x)dx = / (^o S(t) o o To(x)(To(x)dx. 

Then, using To#(To = i^o we apply (|2.5I) to obtain 

/ o S(() o o To(x)(To(x) dx = / o S(t) o $(t)(y)j/o(y)dy, 

JO J A 

since o S(t-) o G i°°(A). Then, since o S(t) G i°°(A) we can use (|6.13p and apply 
(P3|) to get 

/ (^oS(t)oci>(,)(y)z.o(y)dy= / ^ o S(,) (Y)KY) dY. 

J A J A 

Finally, since S(t) satisfies S(t)^i^ = a, we have that 

/ (^oS(t)(Y)KY)dY== / (^(XHX)dX. 
JA Jn 
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Thus, we have shown that 

/ (^(F(t)(x))cro(x)dx = / v3oS(t) o$(t) oTo(x)o-o(x)dx (7.4) 
Jn Jn 

= / </'oS(t) o$(j)(y)iyo(y)dy (7.5) 

J A 

= I ^oS(t)(Y)KY)dY (7.6) 
ip{X)a{X)dX, 



as required. □ 

We now prove that (|5.10l) holds for aU g G [1, oo). 
Proposition 7.4. For any to G [0, r) and any q G [1, oo), 

hm / |F(,)(x) - F(,„)(x)|Vo(x)rfx = 0. 

Proof. By Lemma [7.11 we have that, for any t G [0,r), (|7.ip holds ctq — a.e. in H.. Thus, 
since To#(To — i'q, we see that, for any t,to G [0,t), 

/ |r(t)(x)-F(,„)(x)|Vo(x)dx= / |S(t)0$(,)oTo(x)-S(,„)0$(,„)oTo(x)|Vo(j(W 
Jn Jn 

|S(t) o$(t)(y) - S(t„) o$(t„)(y)|'i/o(y)rfy (7.8) 

|S(t) o$(i)(y) - S(i„) o$(t)(y) + S(f„) o$(t)(y) - S(fo) o (y) | %o(y) dy (7.9) 

/ |S(t)0<D(,)(y)-S(,„)0$(,)(y)|%o(y)rfy (7.10) 
Ja 

+C [ |S(t„) o $(,)(y) - S(t„) o $(t„)(y)|%o(y) dy (7.11) 



=:C(/i+/2). 

Firstly, we show that Ii — as t — !> to- Note that, using (|6.14p and Lemma fe.ll (iv) 
we have that ||i^(t) |lLg(A) ~ II'^o|Il?(a) ^'^^ ^ ^ so that ||i^(t) ||^,(-a) independent of t. 

Let r and r' be conjugate exponents (i.e. - + ^ = 1), with 1 < r < cxd. Then, we can use 
(j6.13[) and Holder's inequality to estimate 

h = [ |S(t)0$(,)(y)-S(t„)0$(,)(y)|%o(y)dy (7.12) 
J A 

= / |S(,)(y)-S(t„)(y)|%(y)dy (7.13) 

J A 

^ |^|S(,)(y)-S(,,)(y)|^'-' {^IKy)r ^y} " (7-14) 

= l|S(t) -S(t„)||J^, (7.15) 
= ||S(t) - S(fo)||^^, lli^ollr as t to- 
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Next, we show that I2 as t ^ Iq. Since S(t) G fl for each t and for v — a.e. y, then, 
by the dominated convergence theorem, it remains to prove that for every Iq, 

S(to)°$(t)(y)-S(t„)o$(i„)(y)^0 ast^to (7.16) 

for v — a.e. y G A. First we note that, since $(t) is measure preserving, then it follows from 
Lemma lOfz). and the fact that w G L°°([0,t) x R^) by jOl), that 

*(t)(y) ^ *(to)(y) ast^to 

in [0,r] for 1/ — a.e. y G A. If y is such a point and if, in addition, $(((,) (y) is a point of 
continuity for ^[to)i then convergence in (j7.16l) holds at y. Since $(to) is measure preserving, 
it follows that $(to)(y) is a point of continuity for S(fQ) for v — a.e. y. Thus, (|7.16p holds 
for J/ — a.e. y G A. □ 

Lemma 7.5. Let Z be defined as in i5.12\) with F defined as in |y. j[ ). Then, for all 
<G [0,t), 

Z(t)(x) = $(,)oTo(x) 

/or (To ~ fl-e- X G r2. 

Proof. Using (|5.12p . we have that Z((-) = T((') o F((-). Therefore, we need to justify the 
following formal computation: 

T(t) o F(t) = T(t) o S(t) o o To = $(t) o To 

since, by p.l6p . T(t) ° ^(t) is the identity on the support of i'. 

Now we make this argument rigorous. Since T exists and is unique a — a.e. in for 
every t G [0,t), we have that T exists and is unique in i7 \ for every t G [0,r), where 

is a Borel subset of fl with cr[iV^] — 0. Then, by Proposition 17.31 we have that 



CTO 



f]nF7j|(Ari)l = / CTo(x)dx= / a{X)dX = 0. 



Now, using Lemma mi we conclude that 

Z(t)(x) =T(t)oS(t)o$(t)oTo(x) 

for xen\N where ctq [N] =0. _ 

Let M = {y G A : T(t, S(<,y)) ^ y}. Then M is a Borel set. Now, the proof of the 
lemma will be completed if we show that 

ffo [{x G 17 \ TV : o To(x) G M}j = 0. (7.17) 

From Theorem 13.51 {Hi) we have that, for any t G [0, r), 

T(t) o S(t)(y) = y for - a.e. y G A. 

Thus, we have 

Jm 
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for any t G [0,t). Therefore, using that (7o[N] = 0, which impUes that To#(To = i^o for all 
X e 57 \ A^, and also using (I6.13|) . we obtain for any t G [0, r) 



{xer!\iV : $(t)oTo(x) gM} - ctq {x G 17 \ TV : To(x) G $^,)(M)} (7.18) 



/ _ cro(x)dx (7.19) 

Jt-1(*-1(M)) 

i^o(y)dy (7.20) 



*-;(M) 

/ i.(y) rfy = 0, 
Jm 



IM 

as required. □ 

We now show existence of the map F* from Definition 15.21 iiii). 
Proposition 7.6. The map F as defined in |7.i[ j satisfies property (iii) of Definition \5.^ 
Proof. As with Lemma [7?T1 we can show that for every t G [0, r) the expression S(o) ° ^"((-1 ° 



-it) 



(x) is defined for a — a.e. x G 17, and the map 



is Borel: 



Since T exists and is unique a — a.e. in fl for every t G [0,t), we have that T 
exists and is unique on £7 \ for every t G [0,t), where iV-'^ is a Borel subset 
of fl with cr[7V^] = 0. Also, since S(o) exists and is unique i^o — a.e. on A, we 
have that S(o) exists and is unique on A \ Nq, where A^g is a Borel subset of A 
with vo[Nq] — 0. Then, we have that S(o) o $*^^ o T(j)(x) is defined for aU 

X G n\{N^UM), 

where 

M={Xe{n\N') : cI,*(i,T(t)(X)) ga^o'}- 

Note that, from its definition, M is a Borel set. We must now show that 
a[M] = for every t G [0,r). Fix t G [0,r). Then, using that T(t)#(T = i/ and 
thus T(t)#(T = for all X G 17 \ A^^, and using (j6.13p . we can apply (|2.5p and 
compute 

a [M] - a [{x G f7 \ A^i : T(,)(x) G $(t)(Aro')}] (7-21) 
- / (T(x)dx= / iyiy)dy (7.22) 



j 

Jn: 



vo{y) dy = 0. 



Thus, we can define F^^^ = S(o) ° ° T{t) and F* is a Borel mapping. 

We can now prove that property (Hi) of Definition 15.21 holds. Since F(t)#cro = cr, we 
have that F*^^ o F(t) (x) = S(o) o ^*^^^ o T(t) o F(t) (x) for cr — a.e. x G 17. Then, using Lemma 
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I7.5[ we get F^^j o F(f)(x) = S(o) o (J)*^^ o <J)(j) o To(x) ctq - a.e. in V,. Since, by Lemma lO 
(it;), o$((j(y) = y for v — a.e. y and thus for i/q — a.e. y £ A, and since To#cro = ^'o, we 
have $*jjO$(t)oTo(x) = To(x) for ao-a.e. x e f^. Thus, F*j)oF(t)(x) = S(o)oTo(x) =x 
for cr — a.e. x G f2 by Lemma l3.5l (Hi). 

By a similar argument, we we have that F(t) o F*^^ = x for cr — a.e. x G fi. □ 

Finally, we show that property (iv) of Definition 15.21 holds for F defined in (|7.ip . 

Proposition 7.7. Let F 6e defined as in |7.i[ ). Then, equality ^5.14^ holds for any Lp G 
CcH(0,t) X f^;R3). Moreowr, we have that Z(-,x) G iyi'°°([0, r)) /or ctq - a-C- x G f^, and 

Proof. From the definition of the Lagrangian flow $ in Lemma 16. f I we have that 

$(i,y) = y + / w(s,$(,)(y))ds 
Jo 

for i/Q — a.e. y G A and every t G [0, r). Thus, this equality holds for all y G A \ TV where 
1^0 [A''] = 0. Since To#tTo = vo it follows that 

ao [n n To i(A^)] = / (To(x) dx = / i.o(y) dj = 0. 

Thus, for (Tq — a.e. x G 51 and every t G [0, r), we have 

To(x)) - To(x) + / w(s, $(,)(To(x))) ds, (7.23) 
Jo 

where we have replaced w(s, (To(x))) by w(s, $(s)(To(x))) based on ([^ . (p?71) . Mul- 
tiplying (I7.23P by cro(x) and by dtfitjX), where (p G C^([0,r) x R^), and then integrating 
we obtain 

/ dMt,x)^{t,To{xj)ao{x) dtdx = I dt(p{t,x)To{x)ao{x)dtdx (7.24) 
J[o,T)xn J[o,t)xsi 

+ / (To(x)9t(^(i,x) /" w(s,$(,)(To(x)))dsdtdx. 
J[o,t)xo Jo 

Now, in the right-hand side, we perform the integration with respect to t in the first integral 
and integrate by parts with respect to t in the second integral to obtain 

/ at(^(t,x)$(t,To(x))(7o(x)dtdx = - / (p(0,x)To(x)cro(x)dtdx (7.25) 
Jlo.T)xn Jn 

(p(t,x)w(i,$(t)(To(x)))cro(x)dt(ix, 



'[0,T)xn 

'[0,r)xn 

(7.26) 

where we have used that (p{t,x) = due to its compact support. Note that, by (|3.15p . 
(17.11) and Lemma [7.51 we have 

w(i,$(4)(To(x))) = /eo.eax [$(,)(To(x)) - S(i)($(,)(To(x)))] = Aoreg x [Z(t, x) - F(t, x)] 

(7.27) 
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for CTo — a.e. x € and every t S [0,t). Substituting (|7.27p into the right-hand side of 
([7^ and using Lcmma[73]to replace To(x)) by Z{t, x) in the left-hand side of ([7?^ . 
we obtain 

/ dt(p{t,x)Z{t,x)ao{x)dtdx = - ip{0,x)To{yi)<7o{x)dtdx (7.28) 
J[o,T)xf2 -/n 

- / ^(t, x) /cores X [Z(t, x) - F(t, x)] CTo(x) dtdx, 

J[0,t)xO 

and rearranging gives (|5.14p . 

Finally, Z(-,x) e W^^'°°([0, r)) for ctq — a.e. x G il follows from Lemma [7751 and Lemma 
KT\(i). Then, (fT^ . (i7:?7l) and Lemma 1731 implv (|OT|) . □ 

Now the properties of (T,cr) in Theorem 14.11 and the properties of F proved in Propo- 
sitions [721 EISl El El EZl imply Theorem [531 



8 Conclusion 

The main result of this paper is the proof of existence of weak Lagrangian solutions of the 
fully compressible semi-geostrophic equations with rigid boundary conditions, in the original 
formulation with variables expressing physically relevant quantities. This result is stated 
in Theorem 15.51 and can be considered as the conclusion of the analysis of the problem 
given in [S], where an existence result is proved but only in the so-called dual formulation. 
We have also proved that, if additional regularity of the flow could be assumed, this weak 
Lagrangian solution would determine a weak (Eulerian) solution of these equations. 

In addition to the main result, we have given an alternative proof, based on recent 
results of Ambrosio and Gangbo on Hamiltonian ODEs in spaces of probability measures, 
of the previous result on the existence of weak solutions of the dual formulation of the 
equations. 
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